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Background 



One of the most 
distinctive signals of the 
Madden-Julian Oscillation 
(MJO) is the upscale 
development and 
organization of convection 
in the Indian Ocean. 

Dynamics of the MJO 
(DYNAMO) campaign 
occurred in late 2011 - 
early 2012 to investigate 
this genesis stage. One of 
the best non-satellite wind 
datasets ever obtained 
over the ocean. 


The Cyclone Global Navigation Satellite System (CYGNSS) mission can exploit 
this dataset to better understand the performance of the satellite constellation in 
regions of deep convection, in particular for characterizing the MJO. 




Three Scientific Objectives 


1 . Produce a high-resolution surface wind dataset for 
multiple MJO onsets using WRF-assimilated winds and 
other data from DYNAMO. 

2. Use the DYNAMO datasets, along with available 
scatterometer observations, to study the causes and 
impacts of wind variability at spatial and temporal 
scales finer than those planned to be provided by 
CYGNSS, and the implications of these processes for 
CYGNSS observations. 

3. Using a simulated CYGNSS dataset for the MJO, 
perform observing system simulation experiments to 
determine the benefits of CYGNSS for improving 
scientific understanding and forecasting of the MJO, 
particularly its genesis over the Indian Ocean. 



Objective #1 
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Objective #1 


WRF Model Setup 

• Advanced Research WRF 

• A: 9-km resolution Indian Ocean domain 

• B: 3-km DYNAMO quadrilateral domain 

• C: 1-km high-resolution domain focused on Revelle 

• 40 sigma levels (more levels in lower troposphere) 

• Nested 9-3-1 km runs, plus separate 3- and 1-km runs 

• Background runs successfully resolve mesoscale features for MJO events 

• Focus on late October, late November, late December 2011 MJO onsets 



Objective #1 


WRF 3D VAR Progress 

• Finalizing WRF runs for B/R matrices generation 

• Assimilation tests underway 

Assimilation Datasets 

• Radars - TOGA, Mirai, S-PolKa 

• Soundings 

• Surface, Ship, and Buoy observations 

• ASCAT/OSCAT 

• HRDL 
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Objective #1 

Preliminary Radar Data Assimilation Results 

2011-11-24 0800 UTC 



• Results suggest large 
impacts from assimilating 
just TOGA data 

• Improvements in precip 
system structure seen 
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NOAA High-Resolution Doppler Lidar (HRDL) 


Objective #2 


On Revelle for Cruises 1-3 (1 September - 6 December 
2011) 

HRDL scanning ability provided 20-min averaged vertical 
profiles of wind speed and direction from 12.5 m to ~2000 m 
12.5-m gate used to compare HRDL to scatterometers 
(OSCAT and ASCAT) 

Examine relationship between winds and mean square slope 
as Richardson number varies 


HRDL RV Revelle DYNAMO 2011 - Wind Speed (m/s) and Direction Profiles 
00:00 10/30/11 to 12:00 10/30/11 UTC 
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ASCAT/HRDL 


OSCAT/HRDL 






Preliminary Results 

• Richardson number always negative for DYNAMO (unstable) 

• ASCAT matches HRDL better than OSCAT 


CYGNSS E2ES Test Run 


WRF simulation of November 
2011 MJO onset 

11/25/2011 Background 10-m 
Winds 

View of soon-to-be JWTC 
Tropical Storm 05A 


24-h, 1-min resolution 



Longitude (deg E) 


CYGNSS Wind Speed (m/s) 


** * 


Objective #2 


Input Winds and Simulated CYGNSS Data 
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CYGNSS E2ES Python Toolkit 


Objective #2 


• Supports reading/analyzing/plotting of E2ES input and output data 

• Complex plots and analyses in just a few lines of code 

• Can isolate single satellites or specific time periods 

• Object-oriented, scalable 

• Model for Python-based analysis of future CYGNSS data? 


In [17] s 


In [18]: merged_data = pygnss . CygnssMultiSat ( output_L2_f iles ) 

display = pygnss .CygnssL2WindDisplay(merged_data) 
display . histogram_plot ( bins=30 , axis_label_f lag=True ) 



lgist = lighten_cmap(plt.cm.get_cmap( 1 gist_ncar ' ) ) 
fig = plt.figure(figsize=(8,8) ) 

input_display = pygnss . InputWindDisplay ( input_f iles [ 11 ] ) 
m = input_display .basemap_plot(title= ' Satellite #3 1 , cmap=lgist) 
ss_display = pygnss .CygnssL2WindDisplay(merged_data. satellites [ 2 ] ) 
ss_display.specular_plot(cmap='gist_ncar' , colorbar_f lag=False, edge_f lag=True, 
title_f lag=False / basemap=m) 


Satellite #3 
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Summary and Future Work 


• DYNAMO dataset nearly prepped for assimilation into WRF. 

• WRF domains established and preparatory 3DVAR assimilation 
matrices nearly finalized. 

• Initial assimilation testing suggests realistic background runs receive 
large impacts from radar data. 

• HRDL/scatterometer comparison suggests good overall performance, 
with ASCAT (C-band) better than OSCAT (Ku-band). Continue 
examining Richardson number relationship to MSS. 

• E2ES working on WRF output, Python module developed. 

• Objective #3 (OSSEs focused on MJO) - Work planned latest this 
year into next. MJO Nature Run? 



